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Abstract. Multi-connected universe models with space 
identification scales smaller than the size of the observable 
universe produce topological images of cosmic sources. 
We generalise to locally hyperbolic spaces the crystallo- 
graphic method, aimed to detect the topology from three- 
dimensional catalogs of cosmological objects. Our new 
method is based on the construction of a collecting-cor- 
related-pair technique which enhances the topological sig- 
nature and can make it detectable. The main idea is that 
in multi-connected universes, equal distances occur more 
often than by chance. We present an idealised version of 
this method as well as numerical simulations, we discuss 
the statistical relevance of the expected signature and we 
show how the extraction of a topological signal may also 



and t he pr oceedings of the workshop Cosmology and topol- 
ogy ( |l998|) for the latest developments) many methods 
using either the cosmic microwave background (see e.g. 



Stevens et al fll993|), Co rnish et al. (|l998a|; |l998b|; |l998c| ) 
Levin et al. ( |l997j ; |199S| ) and Uzan et al. (| 



L998a 



1998b)), 



or discrete sourc es suc h as clusters, quasars, etc. ( see e.g . 
Lehoucq et al. ( 1996| ) and Roukema and Edge ( 1997 )) 
were investigated. The latter class of methods is based 
on the fact that an observer in a multi-connected universe 
is looking at a finite part of its universal covering space 
(Ellis, 1971). We thus expect to see multiple images of 
any object ( Lachicze-Rey and Luminet, 1995| ). This prop- 
erty was first applied to specific objects to provide bounds 



on the size of the universe using our Galaxy (Sokolov 



lead tc a precise determination of the cosmological pa- 1971; Fagundes and Wichoski, 1987), the Coma cluster 



rameters. A more realistic version of the method which (Gott, 1980; Roukema, 1996) and quasars (Roukema and 



takes Eccount of uncertainties in the position and veloc- 



ity data is then discussed. We apply our technique to 
presently available data, namely a quasar catalog, and 
discuss the significance of the result. We show how the 
improved crystallographic method both competes with 
the two-dimensional methods based on cosmic microwave 
background analyses, and suffers from the same drawback: 
the shortage of present observational data. 
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1. Introduction 

Recently there have been many advances in the develop- 
ment of methods to detect or constrain the topology of 
the universe (i.e. of its spatial sections). Since the revival 
of cosmic topology (see Lachieze-Rey and Luminet (1995) 



Edge, 1997). R was then applied statistically through the 
so-called crystallographic method (Lehoucq et al., 1996), 



the applicability of which was discussed by Fagundes and 



Gausmann (1997) in Euclidean manifolds and by Lehoucq 
et al. (1999) in hyperbolic and elliptic manifolds, see also 
Gomero et al. (1998), Fagundes and Gausmann (1998). 
The crystallographic method is a statistical one which was 
framed as to extract characteristic lengths associated with 
the topology of the universe, and thus which gets around 
the problem of recognizing object s at different ages and 
orientations (Lehoucq et al., 1996). 



Today, multi-connected universes with locally Eu- 
clidean spatial sections are well constrained both by the 
absence of a cut-off in the angular power spectrum of the 
cosmic microwave background temperature anisotropies 
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(Stevens et al., 1993), and by the crystallographic method 
( Lehoucq ct al., 1996| ; Uzan et al., 199G) which puts a 
lower bound on the characteristic size, L say, of Eu- 
clidean space models to L > 3000/i~ 1 Mpc (with h = 
iJo/lOOkms^Mpc- 1 , H Q being the Hubble constant). 
But there is still no reliable constraint on a universe with 
locally hyperbolic spatial sections either from the cosmic 
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microwave background methods (Inoue, 199S; Bond et al 
1998|)j>r from the clusters and quasars catalogs (Lchoucq 
et al. , 1 99E ) . Thus , the quest for cosmic topology will prob- 
ably be longer than expected. 

The goal of this article is to improve the crystallo- 
graphic method so as to make it suitable for the detection 
of topology in locally hyperbolic spaces. For that purpose, 
we first review the crystallographic method in Sect. ^], 
comment on its applicability and also on recent attempts 
at its improvement and generalisation. This leads us to 
propose a new method to detect the topology using a cat- 
alog of discrete sources (Sect. |), aimed to enhance the 
topological signal by collecting all correlated pairs of im- 
ages. We first apply the technique to simulated catalogs 
(Sect. ||) in order to study its statistical relevance, then 
to a real catalog of more than 11,000 quasars (Sect. |J). 
We discuss this result after having studied all the uncer- 
tainties in the implementation of the method and in real 
data. 

2. Ups and downs of the crystallographic method 

Assuming that the four-dimensional manifold M. describ- 
ing the universe is globally hyperbolic ( Hawking and Ellis 



1973), which is the case since we reduce our study to uni- 



verses which can locally be described by a Friedmann- 
Lemaitre spacetime, it can be split as M. = R x E 
where the spatial sections E are locally homogeneous and 
isotropic. 

Now, if the universe is multi-connected, E can be writ- 
ten as 



E = X/T, 



where X is either the Euclidean space {E 3 ), the spheri- 
cal space (S 3 ) or the hyperbolic space (H 3 ), and T is the 
holonomy group, namely a finite discrete subgroup with- 
out fixed point of the full isometry group of X (Lachieze 
Rey and Luminct, 1995| ) . The determination of the uni- 
versal covering X will be given by the knowledge of the 
cosmological parameters such as the density parameter fio 
and the cosmological constant A, which will help us to de- 
termine the local geometry (since Einstein's equations are 
partial differential equations, there is locally no difference 
between a multi- and a simply-connected universe). 

The three-dimensional manifold E can be conveniently 
described by its fundamental domain, which is a polyhe- 
dron with 2K faces identified by pairs through the ele- 



ments g of the holonomy group T (Lachieze- Rey and Lu- 
minct, 1995| ). A classification of these manifolds can be 



found in Lachieze- Rey and Luminet (1995) for Euclidean 
and elliptic manifolds, and by using the software SnapPea 



(tittp: //www. northnet.org/weeks) for hyperbolic man- 



ifolds, see also Thurston (1979) and Wolf (1984) 



The crystallographic method (Lchoucq et al., 1996) 



relies on a property of multi-connected universes accord- 
ing to which each topological image of a given object is 



linked to the other one by the holonomies of space which 
are indeed unknown as far as the topology is not deter- 
mined. The only thing we know is that these holonomies 
are isometries. For instance in a 3-torus universe, to each 
holonomy g is associated a distance A g , equal to the length 
of the translation by which the fundamental domain is 
moved to produce the tessellation in the covering space. 
Assuming the fundamental domain contains N objects 
(e.g. galaxy clusters), if we calculate the mutual 3D-sepa- 
rations between every pair of topological images (inside 
the particle horizon), the distances X g will occur at least 
N times for each copy of the fundamental domain, and all 
other distances will be spread in a smooth way between 
zero and two times the horizon distance. In a pair sepa- 
ration histogram which plots the number of pairs versus 
their 3D-separations, the repetition scales X g should thus 
produce spikes. Simulations in locally Euclidean manifolds 
showed that the pairs between two topological images of 
the same object drastically emerge from ordinary pairs 
(Lchoucq et al., 1996) in the histogram. 



As it was recently explained ( Lchoucq et al., 1999 ), 
two kinds of pairs can create a spike, namely (we keep 
the notation and vocabulary introduced in Lehoucq et al. 
( |1999D ) 



1. Type I pairs of the form {g(x), g(y)}, since 

dist[g (x), g(y)] — dist[x,y] for all points and all ele- 
ments g of r. 

2. Type II pairs of the form {x, g(x)} if dist[x, g(x)\ = 
dist[y, g(y)] for at least some points and elements g of 



(1) It has been shown (Lchoucq et al., 1999; Gomero et al 



1998) that in hyperbolic manifolds, type II pairs will not 



exist since the holomies are not Clifford translations (i.e. 
they are such that dist[x, g(x)} depends on x). 

Indeed, type I pairs are always present whatever the 
topology but their number is roughly equal to the num- 
ber of copies of the fundamental polyhedron within the 
catalog's limit. This number is too small to create a spike 



in the pair separation histogram (Lehoucq et al., 1999) 
and depends on the cosmological parameters since a cata- 
log of redshift depth z will contain all the objects located 
within the 2-sphere centered on the observer and of ra- 
dius xi z )i the radial distance obtained by integration of 
the null geodesic equation up to a redshift z. 

Two possible solutions were proposed to improve the 
crystallographic method in order to detect the topology 
and thus to extract the signature of type I pairs in cata- 
logs. 



1. 



Fagundes and Gaussmann (1998) developed a variant 
of the crystallographic method where they subtract the 
pair separation histogram for a simulated catalog in a 
simply-connected universe from the observed pair sep- 
aration histogram (with the same number of objects 
and the same cosmological parameters). The result is 
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"a plot with much oscillation on small scale, modu- 
lated by a long wavelength quasi-sinusoidal pattern". 
The statistical relevance of this signature still has to 
be investigated. 
2. Gomero et al. ( [1998 ) proposed splitting the catalog 
into "smaller" catalogs and averaging the pair separa- 
tion histograms built from each sub-catalog to reduce 
the statistical noise and extract a topological signal 
from what they call "the non-translational isometries" . 
The feasability of this method has not yet been demon- 
strated. 

We now develop a new method based on the property 
that whatever the topology, type I pairs wi ll always exist 
Indeed, as we have shown in Lehoucq et al. (Lehoucq et al. 



199S ) , this does not lead to any observable spike in the pair 



separation histogram. We thus have to improve the crys- 
tallographic method to enhance the topological signature 
by "collecting" all these correlated pairs together. This can 
be achieved by defining what we call a collecting correlated 
pair method (hereafter CCP-method). Such an approach 
(described in detail in the next section) will not be able 
to determine the exact topology but will give a signature 
of the existence of at least one compact spatial dimension 
on sub-horizon scale, which is indeed a first step toward 
the determination of the topological scale of the universe. 
Once such a signature is obtained, many routes are possi- 
ble to determine the exact topology. We shall briefly refer 
to them in the final discussion. 

3. How can we use a 3D-catalog to detect the 
topology of the universe ? 

3.1. Basic idea 

As stressed in the former paragraph, type I pairs will al- 
ways exist in multi-connected spaces as soon as one of the 
characteristic lengths of the fundamental domain, such as 
the injectivity radius rj n j (see, e.g., fig. 10 in Luminet and 



Roukcma ( 1999] )), is smaller than the Hubble scale. Defin- 
ing 5i| 1 < i <2K as the 2K generators of T and referring to 
x as the position of the image in the universal covering 
space X, we have: 



1. Vx,yei,V s er, 
dist[g(x),g(y)] = dist[x,y]. 

We will refer to these pairs as xy-pairs. 

2. VxeX,Vgi,g 2 er, 



(2) 



dis%i(x),gi o g 2 (x)} = dist[x, g 2 (x)]. (3) 

We will refer to these pairs as xg(x)-pairs. 

Both the ccy-pairs and the x<7(a;)-pairs are type I pairs. 

To collect all these pairs and enhance the topological 
signal, we define the CCP-index of a catalog containing 
N objects as follows. 



1. We compute all the 3D-distances dist[a;, y] for all pairs 
of points within the catalog's limit. 



2. We order all these distances in a list d. 



l \Ki<P 



, where 



P = N(N — l)/2 is the number of pairs, such that 
di+i > di. 

3. We create a new list of increments defined by 

Vie [1...P-1], Ai = d i+1 -di (4) 

(keeping all the equal distances, if any, in the list). 

4. We then define the CCP-index K as 



where N = card({i, A 4 = 0}), so that < K < 1. 



(5) 



With such a procedure, all type I pairs will contribute 
to A/. For instance, if a given distance appears 4 times 
, it will contribute to 3 counts in J\f. 



Ki<P 



in the list di 

Compared to the former crystallographic method, all the 
correlated pairs are gathered into a single spike, instead 
of being diluted into the noise of the histogram pair sep- 
aration. 

Indeed, in a more realistic situation, one has to take 
into account bins of finite width e and replace TV by 



M e = card({i,Ai € [0,e[» 



(6) 



in the computation of 1Z. The effect of such a binning will 
be discussed in Sect. ||. We now focus on the "idealised" 
version of the procedure and study the magnitude of the 
CCP-index in multi-connected models. 

For that purpose, let us assume that the catalog is ob- 
tained from an initial set of A objects lying in the funda- 
mental domain and that there are B copies of this domain 
within the catalog's limits [B = if the whole observ- 
able universe up to the catalog's limit is included inside 
the fundamental domain). The total number of images is 
N = A(B + 1). Indeed B is not usually an integer but we 
assume it is, in order to estimate the amplitude of 1Z and 
compare it with the result in a simply-connected model. 

We now need to count the number of copies of a given 
pair. 

1. Each xy-pair will be represented B + 1 times, so that 
the A(A-l)/2 pairs contribute to BA{A-l)/2 counts 
in TV. 

2. For x<7(a;)-pairs, let ^i(E, B) be the contribution to M 
when A = 1. For all x £ E, the total contribution of 
the xg(x)-p&irs is then A 2 i/i(E, B). i/i(Y,,B) can be 
computed by considering a catalog generated from an 
initial set of A = 1 object and has the property that 
whatever the manifold S, ^i(S, 0) = and Vi(T,,B) is 
an increasing function of B. 

3. Summing these two type I pair contributions we ob- 
tain: 



AL 



.4 



(A-1)|+^(E,B) 



(7) 
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Fig. 1. Example of the 2-torus with B = 8. We have 
represented the different type I pairs that contribute to 
^i(T 2 ,8). 



Indeed, if the holonomy group T contains Clifford 
translations allowing for type II pairs, or if there are "fake" 
pairs (i.e. such that dist[x, y] = dist[u, v]), Af m i n computed 
from (|7|) will give a lower bound for the true Af. 

The normalised CCP-index (||) follows straightfor- 
wardly. 1Z is a good index for extracting the topological 
signal since 

1. when B = (i.e. when the fundamental domain is 
greater than the catalog's spatial scale), 1Z = 0, 

2. when the number of sources in the fundamental do- 
main becomes important, it behaves as 



n 



+ 2^(E,g) 
(B + l) 2 



.4 



CO. 



(8) 



3. 2. Application to a 2- dimensional flat manifold 



As an illustrating example, we perform the computation 
in 2-dimensional flat manifolds, and more particularly in 
the case of the 2-torus T 2 . Such a space contains Clifford 
translations, however we count only the contribution of 
type I pairs to 1Z. In that case, the only quantity which 
needs to be computed is Vi(T 2 , B). 

Consider the 2-torus T 2 with its eight nearest neigh- 
bours (B = 8) [see Fig. |l| for the notation]. The pairs 

— (x,gi(x)) and (cc, 174(2;)) appear 6 times, 

— (x, gi o g2(x)) and {x, g\ o (74(2;)) appear 4 times, 

- (92og 3 (x), gi og 2 {x)) and (g 2 og 3 (x), 53° 9±{x)) appear 
3 times. 

- (94,(x), 92 °93(x)), (g 4 (x),g 1 og 2 (x)), (g^ix), g 2 o g 3 (x)), 
(93(x),giog 2 (x)) and (g 3 (x), giog 4 (x)) appear 2 times, 

so that ^i(T 2 ,8) = 2x5 + 2x3 + 2x2 + 4x1 = 24. This 
leads to N = 4A(7A - 1) whereas P = 9A(9A - l)/2 so 
that 



— < K < — , VA> 1. 

35 ~ ~ 81' 



(9) 



Using the same technique one could compute the same 
quantities for various 2-dimensional topologies such as hy- 
perbolic n-holes tori. 



3.3. Application to 3-dimensional manifolds 

In the case of 3-dimensional manifolds, the computation of 
z^i(£, B) follows the same lines as in the previous section. 
We just give two examples. 

— 3— torus with its 26 nearest neighbours: We use 

the numerotation of Fig. || and the notation^ 123 = 
gi 92 #3 {%) , = x etc. We find that the pairs 
a: (0,1), (0,2) and (0,3) appear 18 times, 
b: (0,12), (0,15), (0,23), (0,35), (0,16) and (0,13) 

appear 12 times, 
c: (0,123), (0,234), (0,345) and (0,135) appear 8 
times, 

d: (1,4), (2,5) and (3,6) appear 9 times, 

e: (12,45), (15,24), (13,46), (34,16), (23,56) and 
(35, 26) appear 3 times, 

f: (1,45), (1,24), (2,45), (2,15) and 8 other pairs ob- 
tained by permutations appear 6 times, 

g: (6,123), (6,234), (6,345), (6,135), 
(1,456), (2,456), (2,156), (4,156), 
(5,246) 

and (5, 126) appear 4 times, 
h: (24,156), (12,456), (15,246), (45,126), 
(13,456), (23,456), (23,156), (43,156), 
(53, 246) and (53, 126) appear twice, 
i: (123, 456) appears once, 

so that ^i(T 3 ,26) = 3 x 17 + 6 x 11 + 4 x 7 + 8 x 3 + 12 x 
5 + 12 x3 + 12 x1 + 1x0 = 282. The eight families of 
pairs are represented on Fig. |^. In that example about 
80% of the increments are in Af. Again, we recall that 
when A = 1, we only collect type I pairs. 
Weeks manifold with its 18 nearest neighbours : 



(1,245), 
(4,126), 



(13,246), 
(43,126), 



This manifold (Weeks, 1985) is the smallest known 
compact hyperbolic manifold. Its description can be 
found by using the software SnapPea where it has the 
closed census m003(-3,l). Its volume, in units of the 
curvature radius, is 0.94272 and its fundamental poly- 
hedron has 18 faces (see also Appendix A in Lehoucq 
et al. (1999)). It can then be shown after (some tedious 
that 



computations) 
v x (Weeks, 18) = 90, 
which implies that TV = <L4(11A - 1). 



(10) 



The CCP-index 1Z is given as a function of v\, A (the 
number of sources in the fundamental domain) and B (the 
number of copies of the domain within the observable uni- 
verse) by: 



K = 



A[{2vi - B)A - B] 
IVA 2 



(B + l) 2 A 2 -(B + l)A-2' ^ 

In figures ^ and ^ we plot 1Z as a function of the number 
of objects N for the two 3-dimensional examples given 

1 This labelling of the cells is inspired from the labelling of 
the pr incipal axes in (ordinary) crystallography (see e.g. Phan 
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3-torus, B=26 






Fig. 2. Example of the 3-torus with B = 26, convention 
for labelling the copies of the fundamental domain. The 
image of a point in the fundamental domain (at the center) 
is mapped into the cell e.g. 123 by g\ o g 2 o g 3 . 



2 2.5 
log(N) 

Fig. 4. Plot of U(N) for the 3-torus T 3 with B = 26 
copies. 

Weeks, B=18 






2 2.5 
log(N) 

Fig. 5. Plot of 1Z(N) for the Weeks manifold with B = 18 
copies. 



_.„„ . rn „ „ ,.„. r-i- r indicator for the existence of a non-trivial topology on 

Fig. 3. Example of the 3 torus. Ihc different families of , „ 

7 ,, .. Irr ^n . . . . . scales smaller than the catalog s limit, 

pairs that contribute to v\\L , 26) as labelled m the text 

starting from the base point in 456. 

4. Numerical simulations 




above, respectively {T 3 , 26} (N = 27 A) and {Weeks, 18} 
(N = 19 A). 

By comparison, when computing numerically the 
CCP-index in a simply-connected Friedmann-Lemaitre 
model with the same number of objects and copies, one 
gets 1Z(N) ~ with some negligible fluctuations. 

Thus, in this idealised method (which neglects the 
thickness of the bin e), 1Z provides a statistically relevant 



In the previous section we have shown how 1Z is a good 
indicator of the existence of a topological lensing. We now 
test our new method numerically and show how to im- 
plement it. For that purpose, assuming a locally hyper- 
bolic space, we generate a toy cata log of cosmic sources 
as explained in Lehoucq et al. (199E) and we compute the 
CCP-index. 

First, we need to determine the radial distance as a 
function of redshift (since the catalogs provide the angu- 
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lar coordinate and the redshift of the sources), which im- 
plies knowledge of the cosmological parameters. The local 
geometry is given by the Friedmann-Lemaitre hyperbolic 
metric 



ds 2 = -At 2 + a 2 (t) (dx 2 + sinh 2 X dft 2 ) 



(12) 



where a is the scale factor, t the cosmic time and dfl 2 = 
AO 2 + sin 2 Odip 2 the infinitesimal solid angle. 

We compute the relation between the radial coordinate 
distance (x) and the object's redshift z. We start from the 
Einstein equations for the metric (|12|) 



H 2 



K 

~n2 



A 
~3' 



(13) 



p m being the matter density, A the cosmological constant 
and k = 8ttG/c 4 . H is the Hubble constant defined by 
H = a/a with X = dtX. We choose the units such that 
the curvature index is K = — 1. Introducing = A/3H 2 , 
fl m = np m /3H 2 and the redshift z defined by 1+ z = oo/o, 
( |l~3| ) can be rewritten as (see e.g. Peebles ( |1993| )) 

H 2 

—2 - o ro0 (i + z f + n A0 + (1 - fl m0 - «Ao)(l + ^) 2 -(i4) 

tin 



We have used equation ( |13[ ) evaluated today (i.e. at t = to) 
and we have assumed that we were in a matter dominated 
universe, so that p m oc a~ 3 (this hypothesis is very good 
since we restrict ourselves to small redshifts). 

The radius of the observable region at a redshift z is 
given by integration of the radial null geodesic equation 
dx = At /a and reads 



a da 
„ a'a 



lAOX 



(i - a, 



^Ao) 



(15) 



Equation (|l^) is integrated numerically and the result can 
be compared, when £!a = 0, to the analytic expression (see 



e.g. Gradstheyn and Ryzhik (1980)) 
2(1 — Sl m o) 



xO) 



arg cosh 1 



Or. 



The determination of the radial distance requires the 
knowledge of the cosmological parameters f2o and A and 
can be obtained analytically only when A = 0. Thus, if the 
universe is multi-connected on sub-horizon scale, the plot 
of 1Z in terms of f2o and f^A should exhibit a spike only 
when the cosmological parameters have the right value (as 
shown in Fig. [?]). If the cosmological parameters are not 
exactly known, the distance (15) will be ill-determined and 
the topological signature will be destroyed (see Fig. |6|). 
This has two consequences: 



1 






* * 



Fig. 6. When the cosmological parameters are well esti- 
mated, the redshift-distance relationship enables to recon- 
struct the isometries (a). However, when the estimation is 
wrong, all the determination of the radial distances are 
biased and there is no more correlated pairs in the coordi- 
nate space (b). The observer is represented by the square 
sitting at the centre. 



1. One should span the parameters space (f^Oi^A) 
in order to detect the topological signal, plotting 

n(n ,n A ). 

2. If there is any topological signal, the position of the 
spike gives the values of the cosmological parameters 
on the scale of the catalog's limit (see Fig. [?]). 

As a concrete example we proceed as follows. We first 
generate a simulated catalog by choosing the number of 
objects in the fundamental domain (A — 30), the topol- 
ogy (Weeks manifold) and the cosmological parameters 
(e.g. fin = 0.2, J1a = 0.1), and we then use a second code 
to apply the test, drawing TZ in terms of the two cosmolog- 
ical parameters. The result is shown in plot 0. We see that 
the method works pretty well in the sense that there is a 
spike signalling the presence of a topology and determin- 
ing the cosmological parameters. But we also can check 
that a slight deviation in the choice of the cosmological 
parameters makes the spike to disappear. This effect will 
be discussed in the next section. 

Now, if applied with the required accuracy for the cos- 
mological parameters, the absence of signature will give 
the lower bound on the injectivity radius of the universe 
(in physical units) 



SOOO/i-iMpc 



> 



din: 



Qax 2 + (1 - fi - ^a) + 



.(17) 



(16) 

5. Working with real data 



When one wants to apply the CCP-method to real data, 
one faces a number of problems. First, we cannot use a 
zero width bin, one of the reasons being that the sources 
are not exactly comoving. Let us first estimate the preci- 
sion needed for the cosmological parameters when working 
with a bin of width e as defined in equation (^). For that 
purpose, we just assume Q\ — and estimate the preci- 
sion on fio from (|l6|), 



fi 



\/Oq£ + 1 
^tt z + 1 - 1 



e = F(Sh,z)e. 



(18) 
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Fig. 7. Weeks A = 30, n = 0.2, n A = 0.1, z = 3 
(quasars). We check that the topological signal stands out 
only at the right values of the cosmological parameters 
selected for the simulation. 
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Fig. 8. Computation of the CCP-index on a simulated 
catalog of depth z = 3 in a hyperbolic universe model 
with the Weeks topology, fio = 0.3, J7a = 0.1, using a bin 
width e = 10~ 6 . 



Assuming f2 G [0.2, 1[ and z g]0, z max ], it is easy to see 
that at f2o fixed, F(Qq, z) is a decreasing function of z and 
that F(VIq, z) — > oo when z — ► 0. Since F(fio, z max = 3) ~ 
1, we deduce that 



sn 



(19) 



Indeed, using a catalog with a smaller depth z max will 
allow us to use a smaller resolution for the cosmological 
parameters. One has thus to find a compromise between 
depth and resolution. A deeper catalog tests larger topo- 
logical scales, but requires a better accuracy of the cos- 
mological parameters and thus a longer computer time. 
In Fig. ^, we give an example with a bin width e = 10~ 6 . 
The binning produces a "background noise" which was 
absent in Fig. 0. 

As real data we now consider a quasar catalog^ ( Vcron 
Cetty and Veron, 1998j ) containing 11,301 objects up to a 
redshift of z max = 4.897 (but only 20% of the quasars have 
a redshift greater than 3). On figures || we have depicted 
its projection on the celestial sphere and the redshift dis- 
tribution of objects. In Uzan et al. ( |1999| ) the pair sepa- 
ration histogram method, valid only in Euclidean spaces, 
was applied to the same catalog and no topological signal 
was found; this raised the lower bound on the character- 
istic size of Euclidean space to L > 3000 /i _1 Mpc. This 
limit, corresponding to Lq/Rh > 0.5, is of the same or- 
der as the bound Lq/Rh > 0.8 obtained from the CMB 
( fBtevcns et al., 1993| ). 



2 although quasars are not as good standard ca ndels as X- 
ray galaxy clusters, see Luminet and Roukema (1999) for a 
detailed discussion. 



We now apply the CCP-method in the Weeks hyper- 
bolic space model assuming Ha — and fl spanning 
[0.2,0.9]. No topological signal is found either. Does it 
mean that there is no topological lens effect on scales 
smaller than z max ~ 3 ? Not necessarily, because we could 
only apply the test with precisions e = 10~ 7 , 10~ 6 , 10~ 5 , 
and were unable to span the cosmological parameter space 
with the required accuracy given by (|l9|). The computer 
time is one of the main limitations of our technique. It de- 
pends on the number of sources in the catalog. Typically, 
for the quasar catalog, we had to order about 64 millions 
pair increments (^|) and the run took about 3 minutes (on 
a SUN Ultra Enterprise 3000 with a 1 Gbytes RAM) for 
each couple of cosmological parameters (Oqj ^a)- Span- 
ning the full set of reliable values ([0, 1] x [0, 1] with their 
sum smaller than 1) with a given resolution 50,/ fi requires 
a total running time proportional to (5f2/i7)~ 2 . Indeed, fu- 
ture observations will hopefully tighten the range of the 
cosmological parameters real values; as a consequence, we 
will have to span a smaller set of the cosmological parame- 
ters and our simulations will escape the present computer 
time limitations. 

Now, a second limitation comes from the peculiar ve- 
locities of the sources. Assuming a typical peculiar veloc- 
ity v p ~ 500 km s" 1 (see e.g. Dekel ( p.994j )) the comovmg 



position of a given object with respect to the observer is 
shifted by the quantity ot(z) such that: 



a(z) 



r(£) 

> : 

a 



(20) 
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Fig. 9. The quasars catalog: (up) the distribution of ob- 
jects on the celestial sphere in Hammer- Aitoff equal area 
projection; (down) the redshift distribution. 

where t(z) is the look-back time. It can be computed fol- 
lowing the same lines as for the computation of (15), which 
leads to 



t(z) 



"0 



ifc y^AO^ 2 + (1 - QmO - ^A0) + ^ 

It follows that, since x < 1 in the integrals ([l5]) and 



t(z) x{z) 
a c 

so that 

v. 



a(z) < -B. x {z) ~ 10- 3 . X W, 



(22) 



(23) 



where x( z ) is the radial distance of the object given by 
(|l6|). a is the uncertainty in the comoving position of a 
source with respect to the observer. We show below that 
the corresponding uncertainty induced in the pair separa- 
tion distances is much smaller. 



For xy-pairs (see Fig. |H)|-a), the space position x' is 
related to its strictly comoving position g(x) by 



g(x) + ^-^-dist[x,g(x)], 



from which we deduce that 



(24) 



dist[x', y' 



dist[g(x),g(y)] + {xy.V 



v p (x) 



xy 



c distfx, y] 
| v p (x) (dist[x,g(x)] - dis%, g jy)}) ^ 
c dist[x,y] 

Since dist[t/(x), g(y] = dist[x, y] the uncertainties come 
from two factors 

1. the velocity gradient; however, assuming large scale 
velocity flows, two neighbouring images have sim- 
ilar peculiar velocities and the gradient is very 
small. 

2. a term proportional to the difference of separation 
distances, which vanishes if g is a Clifford transla- 
tion . Oth erwise (as seen in Fig. 11 of Lehoucq et 
al. (|1999|) ), for neighbouring points dist[x, g(x)] — 
dist[y, g{y)\ <C r+, where r + is the largest charac- 
teristic size of the manifold [i.e. the radius of the 
smallest geodesic ball containing the fundamental 
domain, see e.g. Fig. 10 in Luminet and Roukema 
( [19991 )1 ■ 

Indeed, we will lose some of the signal from pairs be- 
tween far apart objects. 

For xg(x)-pairs (see Fig. |l^-b) the position of gi(x) is 
shifted by u p (x)dist[x, gi(x)]/c so that 



dist[x 3 ,x' 13 ] ~ dist[g 3 {x),g 3 o gi (x)] 



v p (x) 



xgi(x) 



(dist[x,ffx o g 3 (x)} - dist[x,ffi(x)]) 
dist[g 3 (x),gi o g 3 (x)] 



(26) 



Since dist[g 3 (x), g 3 o gi(x)] = dist[x, gi(x)], it is easy 
to see that 



(21) dist[x 3 , x' 13 ] — dist[x2, x\ 



Vp{x) 



12J 



xg ± (x) x 



(dist[x, g 2 o g 3 (x)) - dist[x, gi o g 3 (x)]) 
dist[x,gi(x)] 



(27) 



This strictly vanishes if, for instance, g 2 — 9i ■ Oth- 
erwise, dist[x,<?2 ° g 3 {x)] — dist[x,gi o g 3 (x)] is much 
smaller than r+ . This can be understood from the fact 
that the look-back times of g 2 ° 53 (x) and g\ o g 3 (x) 
are of the same order and anyway, their difference is 
much smaller than the expected time for a photon to 
wrap around the universe. 

To our knowledge, such a discussion about the uncer- 
tainty introduced by proper velocities on the positions of 
topological images has never been correctly addressed. Al- 
though the uncertainty induced by proper velocities is low, 
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Fig. 10. Due to peculiar velocities, two topological images 
are shifted from their strictly comoving position. We rep- 
resent here the relation between the mathematical images 
[<7i(a;)] and the topological images [x'j] (a) for xy-pairs and 
(b) for xg(x)-pau:s. 



it will cause the CCP-index to be smaller than its theo- 
retical value, since some correlated pairs constructed from 
far apart images will be lost. However, the uncertainty 
could further be reduced if, for instance, we were able to 
obtain a 6-D catalog including the peculiar velocities, in 
order to correct the position from the shift given by (53) , 
assuming that this velocity field has a weak time depen- 
dence. One also has to take into account the uncertainty 
in distance determination and spectroscopic errors (see 
( Roukcma, 1996 ) for a discussion of these effects). 



6. Co nclusions 



In this article we have introduced a new method for de- 
tecting the topology of the spatial sections of the universe 
using 3D-catalogs of discrete sources. The main motiva- 
tion for looking for such a method was the failure of the 
standard cosmic crystallography for detecting the topol- 
ogy of locally hyperbolic spaces. 

Our method is based on the construction of a CCP- 
index as detailed in Sect. ^. The main difference with the 
cosmic crystallography is the collecting process of all cor- 
related pairs in the catalog, which enhances the signal as- 
sociated to the existence of a non trivial topology. We then 
showed the statistical relevance of this new method both 
on analytic computations and on numerical simulations. 

Indeed, all the simulations were based on an idealised 
method applied to idealised catalogs. We then took into 
account more realistic situations, discussing the effect of 
the finite width of the bins, of the peculiar velocities of 
the objects and of their distance uncertainties (Sect. ||). 

To finish, we applied our method to a quasar cata- 
log and found no topological signature. As explained in 
Sect. ||, this does not prove that the topological scale of 
the universe is larger than the catalog's size because we 
were limited by computer time in spanning the full pa- 
rameter space with the required accuracy. 

To achieve this task, we would have to test around 10 6 
couples of cosmological parameters, giving a total com- 
putation time on our SUN workstation greater than five 
years. A reasonable computation time can be obtained by 
increasing the calculation speed by two orders of magni- 
tude at least. Thus, we are now trying to implement our 
method on a Cray T3E parallel computer with 288 DEC 



Alpha processors in order to gain a factor of at least 100 
in computation time. It is worth noting that total scan- 
ning of the parameter space can be done by few hundreds 
computers since couples of cosmological parameters can 
be tested independently. 

Indeed, we have to wait for tighter constraints on the 
cosmological parameters coming from various sides of ob- 
servational cosmology (see e.g. the proceedings of the - 
XXXIII rd Rencontres de Moriond (1998) for a review) in 
order to reduce the parameter space. 

Now, it is clear that our method cannot help us deter- 
mine the exact topology of the spatial sections of the uni- 
verse. It will only provide a signature of the compactness 
of these spatial sections on scales smaller than the catalog 
depth. Indeed, once this information is known, one can 
think of developing a shape recognition code (for triplets, 
quadruplets, etc.) that will enable us to reconstruct the el- 
ements of the holonomy group [note that the group struc- 
ture will allow us to reject false identifications and to infer 
the existence of missing objects in the catalog]. 

This approach can be thought of as the 3-D analog 
of the 2-D circles method using the co s mic microwave 
background data (Cornish et al., 1998a ; Cornish ct al 



1998bt pornish et al., 1998cj ). The same as the COBE data 



does not have the required resolution to exhibit pairs of 
matched circles if the universe is small, and the technique 
also has to wait for future satellite missions such as MAP 
and Planck, the CCP method is presently limited to 3-D 
catalogs with low redshifts. Until now we are limited to 
quasars (z max ~ 3). A recent survey (Chen et al., 1998) 
in the Hubble Deep Field south NICMOS field found 17 
galaxies with redshifts between 5 and 10 and 5 galaxies 
with redshifts above 10, among a total of 323 galaxies. 
Such deep surveys can enable us hope that we may apply 
our method to deeper catalogs in the future. Also, just 
as the circle method will suffer from the degradation in 
the circle match caused by detector noise (see e.g. Fig. 7 
of Cornish ct al. ( ggggg )), our CCP-method suffers from 
degradation caused by velocity and distance uncertainties. 

To finish, let us extrapolate a little. If the spatial sec- 
tions of our universe have any (observationally relevant) 
topological property, the graph of 1Z in terms of the two 
parameters (Oo, CI a) should exhibit a resonance at the 
right value of these parameters (otherwise, the distance 
determination being wrong, no signature could be de- 
tected). Besides the detection of a non trivial topology, 
the CCP- method can help us determine the cosmological 
parameters on the scale of the catalog limit and will, in 
that case, be an interesting tool to use together with the 
other methods (theoretical and observational) designed to 
determine these parameters. Note that the CCP-method 
for determining the cosmological parameters is purely ge- 
ometric, contrary to the ones using, e.g., the cosmic mi- 
crowave background anisotropics which assume a given 
scenario of structure formation. Thus, in the case of multi- 
connected universe, it will help to improve the estimation 
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of these parameters from the CMB. Indeed this is, at the 
time being, only prospective. 

In conclusion, our new crystallographic method is well 
backed up mathematically but suffers from practical im- 
plementation difficulties. It is, however, important to re- 
alise that a 3D-catalog of cosmic sources contains more 
information than previously suspected, and that the joint 
study of topology and cosmology can lead to a better de- 
termination of the curvature parameters. 
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